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Abstract 

The beam splitter and phase shifter, which are the key elements in the 

experiments of light interference, are realized in the motion of trapped ions. 

Some applications, such as the creation of quantum motional states and the 

realization of Mach-Zehnder interferometer, are illustrated. Several detection 

methods of motional states used in the interferometer are also discussed. 
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1. Introduction 

Beam splitter is a very important element in quantum optical experiments which can 
realize the splitting and mixing of the input light fields to produce correlated or entangled 
output states. A number of authors have considered the behavior of the quantum-mechanical 
beam splitter in the past few years [1-6]. In nearly every interferometric apparatus such 
as Mach-Zehnder[l, 7-9], Fabry- Perot [8] and Michelson-type[10] interferometers, the beam 
splitter plays a central role. Both passive and active interferometers [10-12] that have been 
touted to revolutionize ultrasensitive measurements like gravity-wave detection, sensing of 
tiny rotation rates, etc. use beam splitters. In addition, phase shifter, which may be realized 
by changing the length of an optical path or by employing the optical Kerr effect [10, 13], 
is also an important ingredient in most optical interference experiments. 

On the other hand, the developments in the cooling and trapping ions[14] have opened 
a new research for both atomic physics and quantum optics. When ions trapped in an elec- 
tromagnetic harmonic potential, the quantized vibratic motions play a role of boson modes. 
And the interaction between ions and laser fields leads the coupling between the external 
and internal degrees of freedom of the trapped ions[15-16]. As the coupling between the 
vibratic modes and the external environment is extremely weak, dissipative effects, which 
are inevitable from cavity damping in the optical regime, can be significantly suppressed. 
This unique feature thus can be used to realize some cavity QED experiments without using 
an optical cavity. Stimulated by these ideas, many schemes have been proposed for gener- 
ating nonclassical motional states of trapped ions, such as Fock states [16], squeezed state 
[17-19], even and odd coherent states [20-23] and their various kinds of superpositions [24-29]. 
These nonclassical motional states may be used to perform quantum computation [30-32], 
test quantum- mechanical basic principals [32] and study quantum decoherence effects [33]. 
Up to now, motional Fock states, squeezed states, coherent states [34-35] and Schrodinger 
cat [36] for the cm. mode of a single trapped ion have been produced and observed experi- 
mentally. The reconstruction of motional states have also been realized in experiments [37]. 

However, the ability that people coherently manipulate motional states is much lower 
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than that people manipulate light fields. In order to further enrich and develop the con- 
tent of this portion, we here present a method to implement the beam splitter and phase 
shifter in the motion of trapped ions, and then use them to coherently produce some typ- 
ical quantum motional states and to implement the famous Mach-Zehnder interferometer 
in the vibrational motion of trapped ions. Some of the results or ideas may also be used 
directly to implement other vibrational interferential experiments, such as the homodyne 
and Michelson-type interferometer. The construction of this paper is arranged as follows: 
In section 2, we first describe how to construct the two pivotal elements for implementing 
a Mach-Zehnder interferometer — the two-mode beam splitter and phase shifter for the mo- 
tional states of cm. and breathing modes in two trapped ions, where the later is mainly 
based on the dynamical evolution of laser-drived trapped ions first introduced by D'Helon 
and Milburn[38]. In section 3, we then illustrate their applications in the generation of 
quantum motional states. And in section 4, we further use them to implement a Mach- 
Zehnder interferometer for the vibrational motion of trapped ions, where three detection 
ways of motional states are provided. The first two ways are mainly based on the idea of the 
reconstruction of motional states[34, 39], and the third one based on the direct measurement 
of the mean value of a motional operator [40]. A conclusion is given in section 5. 

2. Beam splitter transformation and phase shift in the motion of trapped 
ions 

Beam splitter and phase shifter are two key elements in quantum optical experiments. 
Their realizations in the motion of trapped ions will allow us to implement the similar inter- 
ferential experiments in trapped-ion systems without the use of optical elements. Though 
some papers have already involved this problem [4 1-42], the general constructions of beam 
splitter and phase shifter in trapped-ion systems need further developments. Here we provide 
a method to realize this two important elements in the motion of two trapped ions. 

Following the experimental implementation of reference [43], we consider a two-level 
system which is formed by two 9 Be + ions of mass m trapped in a linear trap. The two 
spectrally resolved ground-state hyperfine levels are defined as the ground and excited states 

3 



of our problem: |F = 2,m F = —2) = \g) and |F = l,m F = — 1) = |e), which form the so- 
called internal degrees of freedom. Assume that the two ions are strongly bounded in 
the y and z directions, but weakly bounded in a harmonic potential in the x direction 
and oscillate around their equilibrium positions: x w = —d/2, x 2 o = d/2. We denote by 
X = (x 2 + Xi)/2, x = (x 2 — Xi)/2 the cm. and breathing mode operators respectively 
which constitute the studied external vibratic motional degrees of freedom. Now we let 
only one of the ions, say the first ion, to be illuminated simultaneously by two classical 
standing wave lasers with the positive frequency parts Ej + ^ = E 0I cos(kiX+(pi)e~ luJlt ~ l ' t ' 1 and 
= Eon cos(kux+(pu)e~ lu ' Ilt ~' 1 ^ 11 ■ Assume that the two lasers have the same amplitudes 
Eqj = Equ = E and effective wavevectors kj = kjj = k, but different frequencies uj and lojj- 
Experimentally, the two standing wave lasers with needed properties can be implemented 
by employing two-photon stimulated-Raman transitions[18, 31]. Of course, in order to let 
a laser illuminate selectively only ion 1 without lighting ion 2, the laser propagation should 
form an appropriate angle with x axis[44] and the projective wave in the direction of x axis 
plays the role of coupling the internal and the wanted external degrees of freedom. Assuming 
that the ion is located at the anti-nodes of both the two standing waves, i.e. ipi = <pn = 0, 
then the Hamiltonian for this system can be written as, 

H — H + Hi n t, (1) 

H = fia + a + ub + b + Uqct z1 /2, (2) 

H int = ^ +1 [(e-^*-^ + e-™"*-*") cosim) + h.c] (3) 

where \x [y) and a (b) are the frequency and annihilation operator of the cm. mode (breath- 
ing mode), uiq the energy difference between ground state \g) and excited state |e) of each 
ion, cr z i and cr + i the usual Pauli operators describing the internal degrees of freedom for ion 
1, Q the Rabi frequency for the ion-laser coupling which is assumed to be equal for the two 
standing wave lasers. 
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We select the two standing waves to excite resonantly both the upper and lower sidebands 
of cm. and breathing modes, i.e. the frequencies of the two lasers satisfying 

LOi = uj - (fx - u), u n = u + (/i - v). (4) 

Under the condition that both cm. and breathing modes are cooled to Lamb-Dicke regime, 
we can expand the Hamiltonian up to the order terms of rj 2 or r/ 2 which represent the 
lowest couplings between internal and external degrees of freedom. Transforming the above 
Hamiltonian to the interaction picture with respect to H , we have 

Hl nt = n V ri r (a + b + ab + )cr xl (5) 

for one kind of laser-phase selection of (f>j = (pu = and 

H int = -i^VVr(a + b - ab + )a xl (6) 



for another kind of laser-phase selection of (J>i = —(/>n = — vr/2, where rj = kJl/Amfx and 



r] r = k\j\j^mv are the Lamb-Dicke parameters corresponding to the cm. and breathing 
modes, a x \ the x component of Pauli operator for ion 1. The exact propagators for the 
above Hamiltonians are 

U\ = exp{—iQr]r] r t(a + b + ab + )a xl } (7) 

for the first kind of laser-phase selection and 

U2 = exp{—Qr]ri r t(a + b — ab + )a x i} (8) 

for the later kind of laser-phase selection. Note that the interactions of U\ and U 2 occur 
only between the internal states of ion 1 and motional states of the system. Now we prepare 
the initial state of the whole system as 

l# = l#i <8> WI2 ® (9) 

with|'0) 1 = (|e) 1 + \g) l )/\/2 and \ip) 2 (which may be a arbitrary pure state) the internal 
states of the two ions and \ip) cr the initial motional state of cm. and breathing modes. By 



performing unitary propagators U\ or U 2 on this initial state respectively, we find that both 
the internal states of the two ions remain unchanged but the motional states experiences a 
beam-splitter transformation: 

BiWW c , r ^- W |4, r . (10) 

or 

B 2 {9) |V) c>r = e~ WJ y |^) c , r (11) 

respectively, where Bi{9) = exp(— i9J x ) and B 2 {9) = exp(—i9J y ) with 9 = 2Qrjrj r t being a 
real angle. The Schwinger angular-momentum operators [45] are defined as 

J x = ( a + b + ab + )/2, 
J y = ( a + b-ab + )/2i, 
J z = (a+a - b + b)/2, 

which satisfy the usual commutation relations for the Lie algebra of SU(2): [J x , J y ] = iJ z , 
[Jy,Jz] — iJx an d [J z , Jx] — iJy The transformations Bi(9) and B 2 (9) represent rotations 
by angle 9 about —x and — y axes in the angular-momentum space respectively, which would 
be distinguished as n/2— and n— coupling beam splitters[l] for cm. and breathing modes. 
When 9 = n / 2, they are referred to as 50/50 beam splitters. 

The phase shift transformation for motional states of one trapped ion may be realized 
by changing the trap frequency for a fixed time [46]. Here we construct the phase shifters 
for cm. and breathing modes of two trapped ions using the dynamical evolution of laser- 
drived trapped ions first introduced by D'Helon and Milburn[38]. For convenience in the 
next sections, we consider an interaction of a classical standing wave with ion 2 in the above 
two-ion system. Assume that the ion is placed at the node of the standing wave and the 
detuning A = u — uj l is set to excite only the cm. motional degrees of freedom but far 
away from the corresponding sideband resonance (|A| ^> /i). In Lamb-Dicke limit and for 
interaction times much greater than the vibrational period of cm. mode, this interaction 
produces the following conditional phase shift [38,42] 
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U 3 = exp[-i X ta + a(a z2 + 1/2)], (12) 

where \ — ^ 2 ^ 2 /(2A) with 77 the Lamb-Dicke parameter of cm. mode and f2 the Rabi 
frequency for the ion-laser coupling. Pauli operator a Z 2 describe the internal degrees of 
freedom for ion 2, and t is the interaction time. The conditional phase shift U 3 occurs only 
between the internal state of ion 2 and the cm. motional state which can be used to realize 
a phase shift transformation on cm. mode. For example, by preparing the initial state of 
the involved system as |<7) 2 <8> \^) c with \g) 2 the ground state of ion 2 and \ip) c the vibrational 
state of cm. mode, then the action of U 3 on this state produces a phase shift on cm. mode: 

P c (y?) |^) c = exp(^a+a) |^> c (13) 

without changing the internal state of ion 2. The phase-shift angle ip = ytjl may be adjusted 
at will by means of altering the interaction time. 

The phase shift on breathing mode may also be realized in a similar way. In the above 
process, if we use the dispersive standing wave to excite non-resonantly the breathing mode 
instead of cm. mode, in the similar conditions one has 

U 4 = exp[-t Xr tb + b(a z2 + 1/2)], (14) 

where Xr = Vr^-r/ '(2A r ) with r] r = 3 _1//4 r] the Lamb-Dicke parameter for breathing mode, 
Q r the Rabi frequency and A r = u — uj l the detuning frequency which satisfies |A r | ^> u 
this time. By preparing the initial state of the involved system as \g) 2 ® \ip) r with \ip) r the 
vibrational state of breathing mode, then, without changing the internal state of ion 2, the 
unitary propagator C/4 produces a phase shift for breathing mode: 

P r (ip r ) = exp(iip r b + b) (15) 

with the phase-shift angle ip r = x r tj1. 

It is worthwhile to point out that the similar two-mode beam splitter transformation of 
eq.(ll) can also be realized in the motion of one ion trapped in a two-dimensional isotropic 
harmonic potential trap[41]. The merit to use two ions instead of one ion is that it can 
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realizes simultaneously the two kinds of beam-splitter transformations of eqs.(10)-(ll) by 
simply adjusting the initial phases of the applied lasers. Further, combining these two kinds 
of beam-splitter transformations with the phase shifts given by eqs.(13) and (15), we can 
construct the most general two-mode beam-splitter transformation [6]. The constructions of 
the above beam splitter and phase shifter on motional states are simple, which employ at 
most two standing waves to interact with one ion. Thus, it may be hopeful to demonstrate 
in current or future experimental conditions. 

3. Creation of quantum motional states using beam splitter and phase shifter 
As they played roles in the quantum optics, the beam-splitter and phase shifter for 
motional states constituted in the last section may also possess underlying applications. In 
this and next sections, we will show this compendiously through concrete examples. In this 
section, let us first describe the applications in producing quantum motional states. The 
first kind of important motional states to be produced are the entangled number states. By 
preparing the initial motional state in eq.(9) as \tp) cr = |l) c <E> |0) r = |l,0) cr with |l) c and 
|0) r the Fock states of cm. and breathing modes respectively, and performing U\ or U 2 on 
it, we obtain the output states for motional states as 

B 1 (9) |1, 0) Cjr = cos(£/2) |1, 0) c , r ~ i sin(0/2) |0, 1) C>J . , (16) 

B 2 (9) |1, 0) c>r = cos(£/2) |1, 0) cr + sin(0/2) |0, 1)^ (17) 

with 9 = 2Qf]7] r t which can be adjusted by altering the interaction time t. When 9 = ir/2, we 
obtain the maximally entangled single-phonon number states (jl,0) cr — i |0, l) cr ) /a/2 and 
(|1> 0)cr + |0' l)cr) />/2- Instead, if we prepare the initial motional state as \ip) cr = |1, l) cr , 
then the output states are 

BiiP) |1, l) c , r = cos 9 |1, l) Cjr - i sin 9 (|2, 0) cr + |0, 2)^) /A (18) 
B 2 (9) |1, l) c r = cos 9 |1, l) Cjr - sin£ (|2,0) c>r - |0,2) c r ) /y/2. (19) 
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For 9 = 7r/2, the maximally entangled two-phonon number states (|2,0) cr ± |0,2) cr ) /y/2 
have been created. These maximally entangled number states play an important role in 
quantum measurement theory. Similar motional states have been produced for one ion 
trapped in a two-dimensional isotropic ion trap [41]. The initial motional number states 
used above have been already prepared experimentally by simply applying a sequence of 
Rabi 7r pulses of laser radiation on the blue sideband, red sideband or carrier [34] . Therefore 
the method for producing above entangled phonon number states is practical. Once the 
above beam splitter has been realized in experiment, then the creation of the discussed 
motional states will become truth. 

The second kind of important motional states to be produced are the entangled coherent 
states. According to the formula for beam splitter transformation 



BtV) \a) c 



acos(^) -/3sin(^ 



c 



asin(^)+/3cos(^)) (20) 



r 



with \ct) and \(3) are motional coherent states, if we let \tp) cr = N±[\a) c ± |— a) J <8> |0) r with 
N± = [2(1 ± exp(— 2|a| 2 )] -1 / 2 the normalized constants, then equation (11) leads to a pair 
of entangled odd and even coherent states: 

5;^ ± = JV±[|5> c |)9) r ±|-5> c |- ) 9) r ] (21) 

with a = acos(9/2) and (3 = asin(0/2). For a 50/50 beam splitter, i.e. 9 = n/2, one has 
a = (3 = aj\pl and eq.(21) reduces to 

\a; a) ± = N ± [\a) c \a) r ± \-a) c \-a) r \. (22) 

The preparing methods for these types of entangled coherent states in the motion of one 
trapped ion have been provided [22-23, 26]. In the limit of a — > oo and (3 — > oo, the 
above entangled coherent states become Bell states which can be employed to complete 
teleportation of continuous variables [4 7]. In addition, in the initial state of eq.(9), if we let 
\ip) 2 = \g) 2 an d \ip) cr = -^±[|a) c ± |— a) J <S> \(3) r (the internal state of ion 1 may be 
an arbitrary superposition of \g) 1 and |e) 1 ), after performing the transformations of eqs.(8) 



and (12) successively, i.e. U^t^^fo) with t 2 = n / (4:Qr)r) r ) and t 3 = 2n/x, we then obtain 
another type of entangled coherent states (unnormalized) for motional states 

| £ _) c | £+ > r ±| £+ ) c |£_) r (23) 

with e± — (P ± a)/ \[2 while the internal states of the two ions are immune. The odd two- 
mode coherent states in eqs. (22) and (23) are maximally entangled states for two-motional 
modes. 

In the generation of the motional entangled coherent states described above, we used 
the odd and even coherent states as the input states which can be prepared easily by the 
interaction of a ion with lasers [22-23, 26-27]. The advantages to produce motional states 
in terms of beam splitter or/and phase shifter are the simplicity and practicability: it can 
produces very complicated nonclassical motional states through simple networks constituted 
by them. 

4. Mach-Zehnder interferometer for the motional states 

The Mach-Zehnder interferometer is a fundamental interferential device which may be 
applied to the detection of phase sensitivity. In general, the input states of the inter- 
ferometer are Bose modes which creation and annihilation operators d + and d satisfy: 
d + \n) = y/rT+T \n + 1) and d\n) = y/n\n — l) with \n) a number state of the consid- 
ered Bose mode. Leibfried, et al.[46] have recently simulated a Mach-Zehnder interferential 
experiment in trapped ion system which used the internal and external degrees of freedom of 
a trapped ion as the two input Bose modes of the interferometer respectively. The harmonic 
oscillator mode of the cm. motion of the ion along one axis in the trap served as one of the 
Bose modes. In the restricted Hilbert space formed by the two ground-state hyperfine levels 
\F = 2,rriF = —2) = |1) and \F = l,m F = —1) = |0) of the ion, the raising and lowering 
operators cr + and o"_ between the levels have the same action as the annihilation and cre- 
ation operators d and d + of a Bose field respectively. Thus the internal spin motion of the 
trapped ion can simulate another Bose mode in a restricted two-dimensional Hilbert space. 
Here we implement a Mach-Zehnder interferometer by fully employing two vibrational modes 
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of trapped ions — the cm. and breathing modes. It need not restrict one of the input states 
of the interferometer within the two-dimensional number Hilbert space and thus can make 
more general interferential experiments. The sketch for this interferometer is depicted in 
Fig.l. It consists of two 50/50 vibrational beam splitters Si, S 2 and a phase shifter tp on 
cm. mode. The output phonon-number distributions for the two modes are detected by Di 
and D 2 . Note that the Mach-Zehnder interferometer here works in a temporal rather than 
a spatial regime. The cm. and breathing vibrational motions take place in the same spatial 
region which don't propagate along separate paths in practice, and the introduction of the 
mirrors in Fig.l is only for the convenience of illustrating which don't exist actually. In the 
following, we describe concretely the implement of each element. Now the initial state of 
the whole system is prepared as 

l^o) = |^o)i ® \g) 2 ® \in) (24) 

with 1^0)1 — (l e )i — lfl , )i)/v / 2 the internal state of ion 1 and \in) the input motional state 
of cm. and breathing modes entering the interferometer. At the action of eq.(7) for an 
interaction time T = ir / (4Q,r)r) r ) , the evolved state is 

\*h) = |^o)i ® |ff) 2 ® ^' 2)Jx H (25) 

In this way, we realize the transformation of motional state for the beam splitter Si. By 
performing the unitary transformation of eqs.(12) on this state, we then complete the phase 
shift ip on cm. mode and the resulted state becomes 

|^2> = |^o)i ® \g) 2 ® e^e-WW' \in) (26) 

with ip = i] 2 Q 2 t/(4A). Performing again the propagator of eq.(7) for the same interaction 
time T, the final transformation of beam splitter S 2 is completed and the output state 
becomes as 

1^3) = \^o)i ® \g) 2 ® \out) , (27) 

with 
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\out) = e ^mJ* e ^a+a e i(Tr/2)J x |^ ^ 

the output two-mode motional state. According to the idea of Mach-Zehnder interferometer, 
we should measure the phonon-number difference in the two modes at the output of the 
interferometer, or equivalently, the mean value: (J z ) = (out\ J z \out). Assume that the 
input motional state is 

|in) = |0) c ®|a) r , (29) 

i.e. the cm. mode be in a vacuum state |0) c and the breathing mode in a coherent state 
| a) r (usually with large amplitude \a\), then we have 

( J z) = ^ n cos tp (30) 

with n = \a\ 2 the initial phonon number of breathing mode. It is shown obviously that if 
the mean value of (J z ) can be measured by experiment, then the shifted phase tp may be 
determined. The phase uncertainty for this detection may be estimated. From eqs.(28)-(29), 
we have 

( jfj = (out | j\ | out) 

= ^n(l + n cos 2 tp), (31) 

and hence 

(AJ z f = (Jl) - (J z ) 2 = -n. (32) 

The phase error is thus given by 

A AJ Z = 1 

\d(J z )/dtp\ v/nsin^' 

When tp = n/2, the phase error is limited to the minimum \ j\fri which is reverse to the 
square root of the input phonon number. Of course, the sensitivity of an interferometer can 
be improved by designing the input states appropriately [1]. 
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The function of the detectors D± and D 2 is to measure the phonon number difference 
of the output two motional modes, i.e. (J z ). In trapped ion system, the measurement of 
motional states is usually be realized by first mapping these external degrees of freedom onto 
the internals and then making a measurement on the internal states. The measurement of 
internal states may be finished by the method of electron shelving which has the detection 
efficiency nearly close to one [48]. Below, we present three methods to measure (J 2 ). The first 
way is to measure the phonon number distributions of the two output modes independently, 
which may be performed in principle by use of the dynamics of Jaynes-Cummings model 
[34, 39]. A more general measurements of motional states including the determination of 
non-diagonal elements of a density matrix may be realized by the method of quantum-state 
reconstruction [37]. In the above two-ion system, if we let ion 2 interact with a standing 
wave with the ion located in the node of the standing wave, we can then get the Hamiltonian 
as 

H = LU a z2 /2 + /ia + a + isb + b + H int , (34) 

where H int = (Q/2)a +2 exp(—iu L t) sin(A;ir 2 ) + h.c. with k (uj l ) the effective wavevector (fre- 
quency) of the standing wave and x 2 (0+2) the position operator (internal raising operator) 
of ion 2. The meanings of the other symbols are defined as before. By adjusting the laser 
frequency to be resonant with the first red sideband of the cm. mode, i.e. ujl = ujq — in 
Lamb-Dicke limit and under rotating wave approximation, the interaction Hamiltonian in 
the interaction picture can be reduced as the standard JCM 

Hint = A(o- +2 a + cr~ 2 a + ). (35) 

with A = i]Q/2 the coupling strength. For convenience of narration, we expand the output 
motional state of eq.(28) as 

\out) = c mn \m) c \n) r (36) 

m,n 

with \m) c , \n) r the motional Fock states of cm. and breathing modes respectively and c mn 
the complex numbers. In general, \out) is a two- mode entangled state. Performing the JCM 
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interaction of eq.(35) on the interferometer output state of eq.(27) for a interaction time r, 
followed by a measurement of the internal state of ion 2, then the probability that the ion 
2 in the ground state can be expressed as 

P 9 (r) = i[l + 5> m cos(2Arv/^)] (37) 

1 m 

with p m = J2n \ c mn\ 2 the phonon number distribution of cm. mode. A Fourier transforma- 
tion of the measured signal P g {r) allows one to determine the distribution probability p m of 
the output cm. mode[34, 39]. The phonon-number distribution probability p n = J2 m \ c mn\ 2 
for breathing mode can be obtained in a similar way. After obtained the phonon-number 
distribution probabilities p m and p n of the two output modes, we can then easily determine 
the mean value of J z as 

( J z) = \(^Zmp m - ^2 n Pn)- (38) 

m n 

The second way to determine (J z ) is to measure the expansion coefficient c mn , or equiv- 
alently the joint distribution probability p mn = \c mn \ 2 . This can be accomplished in terms 
of the dynamics of two-modes JCM. In the above laser-ion interaction, if we let the ion 2 
locate at the anti-node instead of node of the standing wave, and adjust the frequency of 
the standing wave to excite simultaneously the both lower sidebands of the two motional 
modes, i.e. ojl = ujq — (// + u), in Lamb-Dicke limit and under rotating wave approximation, 
the interaction Hamiltonian in interaction picture is 

H' mt = g{aba +2 + a + b + a- 2 ) (39) 

with g = Qrjrjr the coupling strength. This kind of two-mode JCM interaction has been 
studied extensively in quantum optics [49]. Now if we let eq.(27) as the initial state of this 
JCM, then the dynamical evolution, after inserting eq.(36), leads to 

I^W) = c mn[ c os(gtV mn ) \ m i n i 9) — i sin(gty/mn) \m — 1, n — 1, e)], 

m,n 

where the first, second and third marks in Dirac symbols denote the states of cm. mode, 
breathing mode and the internal state of ion 2 respectively. The unchanged internal state of 
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ion 1 is eliminated. Now, if we again measure the internal state of ion 2, then the measured 
probability in ground state is 

P g (t) = ^[l + Y,PmnCos(2gtV^)]- (40) 

m,n 

In a similar way, we can determine the joint probability distribution p mn of the motional 
state \m,n) cr by Fourier transforming the measured signal P g (t) and then obtain the mean 
value of J z in terms of eq.(38) as well as p m = J2 n Pmn and p n = J2 m Pmn- 

The third way, which is also the simplest way, to determine (J z ) is the direct measurement 
of the mean value of a motional operator[40]. In order to do this, we first perform a 7r/2- 
pulse carrier transition [48] on ion 2 to reset the output state of the interferometer of eq.(27) 
as 

1^4) = |-0o)i<E> |-0o> 2 (8) |out) (41) 

with 1^0)2 = (\d)2 + * l e )2)/v^2 the internal state of ion 2. Then we execute the unitary 
transformation C/3 of eq.(12) on this state. At this moment, the mean value of (a x2 ) in the 
final state is 

(0x2) = (^4 1 U^a x2 U 3 ^4) 

= - (out\ sm(2xta + a) \out) . (42) 

For small interaction time and not very large number of phonons, we can linearize the 
operator as[40] sin(2xta + a) ~ 2xta + a, then the above equation reduces to 

(cx2) = — ^Xt (out\ a + a \out) . (43) 

In trapped ion system, the measurement of (<J X 2) is very simple which may be completed 
by applying a single qubit rotation to transform the internal state of ion 2 from <r z -picture 
to o^-picture, followed by a detection of the ion's fluorescence [48]. Therefore the mean 
number of phonons of the cm. mode in the output state of the interferometer (a + a) can be 
evaluated easily from eq.(43). The mean number of phonons of breathing mode {b + b) can 
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also be evaluated in a similar way. Thus the mean value of (J z ) = [ (a + a) — (b + b)]/2 could 
be obtained. 

5. Conclusion 

In conclusion, we have proposed a method to realize the transformations of beam split- 
ter and phase shifter in the motion of trapped ions and discussed their applications in the 
generation of quantum motional states and in the Mach-Zehnder interferometer. Several de- 
tection methods for the output motional state of the interferometer have also been discussed. 
Further, our research works may be directly applied to implement some other interferential 
experiments of motional states, such as the Homodyne detection and Michelson-type in- 
terferometer [10]. Though the implementation process involves some complicated ion-laser 
interactions, some pivotal techniques such as the basic manipulations of ion-laser interac- 
tion [34-36,43,48] and the reconstruction of motional states[37] have been already realized 
experimentally. We expect that, with the increasing of the experimental precision, our pro- 
posals will become truth and produce an important effects in the future coherent control of 
quantum motional states of trapped ions. 
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Fig.l. Principle of the Mach-Zehnder interferometer. The cm. and breathing modes a 
and b of the two trapped ions are mixed on the first beam splitter Si. After a phase 
shift (p on the cm. mode, they converge at the second beam splitter S 2 . Finally the two 
output motional states are measured by detectors Di and D 2 . 
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